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1. Introduction and summary 

Strongly coupled far-from-equilibrium systems appear in many areas of physics and their the- 
oretical study is notoriously difficult. The gauge/gravity duality is a novel tool to study certain 
classes of strongly coupled field theories, and its extension to far-from-equilibrium situations is 
a very active area of current research. In order to extract useful field theory information from 
computations in the dual gravity theory, it is important to identify field theory observables that are 
meaningful in far-from-equilibrium states, as well as prescriptions for how to compute them in the 
gravity dual. 

In the recent paper [1], to which we refer for a more elaborate introduction as well as references 
to the literature, our collaborators and we studied a notion of time-dependent spectral functions and 
occupation numbers (based on Wigner transforms of two-point functions) in a simple holographic 
model (see [2, 3] for related work in other holographic models). Starting with the vacuum state 
of a strongly coupled 2d conformal field theory, the sudden, homogeneous injection of energy and 
its subsequent thermalization was modeled by an AdS3-Vaidya geometry that interpolates between 
pure AdS at early times and a black brane geometry at late times. In order to compute two-point 
functions of high-dimension operators, we used a geodesic approximation [4]. For timelike sepa- 
rations of the operators, it was necessary to use either a non-standard Euclidean continuation of the 
bulk spacetime, or to use complex geodesies, and both procedures gave the same result. We also 
worked out an approach that went beyond the geodesic approximation, but will not discuss it here. 

In the present contribution, we extend some of the results of [1] in two ways. First, instead 
of stalling with the vacuum, we start with an initial thermal state. After a sudden, homogeneous 
injection of energy, the system will evolve towards a thermal state with higher temperature. Second, 
as in [2], we extract from our results a time-dependent effective temperature, which evolves in a 
non-monotonic way from the initial to the final temperature. 

In section 2, we review the notions of time-dependent spectral function [3, 1] and time- 
dependent effective temperature [2]. In section 3, we present the AdS-Vaidya model that describes 
the evolution from a black brane with an initial temperature to a black brane with a larger final 
temperature via the injection of a shell of null dust. In this setting, section 4 describes the com- 
putation of retarded and time-ordered two-point functions using a geodesic approximation and a 
non-standard continuation to Euclidean signature. In section 5, these two-point functions are used 
to extract time-dependent spectral functions and effective temperatures, which are shown to inter- 
polate between the initial and final equilibrium values. 

2. A notion of time-dependent spectral function and effective temperature 

We probe the thermalisation of the field theory with time-ordered (Feynman) two-point func- 
tions, 

iG F (t 2 ,x 2 ;h,xi) = (T{0(t 2 ,x 2 )0(t\,x\)}) 

= e(t2-h)(0(t2,X2)0(h,X l )) + G(t l -t2)(0(tl,X 1 )0(t2,X 2 )), (2.1) 

and retarded two-point functions (which are real in position space) 

iG R (h,X2;ti,xi) = O(t 2 -ti){[0(t2,X2) t 0(h,x 1 )]). (2.2) 
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From these definitions, using the hermiticity of the operator &(t,x), we can deduce the relation 

G R (t 2 ,x 2 ;ti,xi) = d(t2-h) [G F (t 2 ,x 2 ;ti,xi) + (G F (t 2 ,x 2 ;ti,xi))*] , (2.3) 

which we will employ later on. When the field theory is in equilibrium, the system is translation 
invariant in space and time. The correlators are then functions of the differences only, i.e. 

G R . F (t 2 ,x 2 ;tuxi) = G RiF (t 2 -h,X2-xi;0,0) = G R . F (t,x), (2.4) 

with t = t 2 — t\ and x = x 2 —x\, and their Fourier transforms are given by 

G R)F {(Q,k)= dt dxe im e~ ikx G R . F (t,x). (2.5) 



In momentum space, the retarded correlation function is in general complex valued, and its imagi- 
nary part defines the spectral function, 

p{(0,k) = -21m G R (w, k) . (2.6) 

When the system is at thermal equilibrium, the retarded and time-ordered correlation functions are 
related by the Fluctuation-Dissipation Theorem, 

[\+2n(co)]p(co,k) = -21mG F ((0,k) , (2.7) 

where the occupation number equals the Bose-Einstein distribution: n(a>) = (e m l e — l) -1 for some 
temperature 6. Note that using the previous expressions, we have 

taG , ( .) (—g±l) , (2.8) 

v 7 co=0 

where 

/+<*> f+<>° 

— G R , F {(0,k)= dte' at G RjF {t,x = 0) (2.9) 
oo J — co 

is the momentum average of the correlators. 

In a thermalising field theory, time translational invariance is generically broken. Therefore as 
in [1] we introduce an average time T and a relative time t by 



J = ^ \ tl =T- 

t=t 2 -h 1*2 = T + j 



and we define time-dependent coiTelators by keeping T fixed and Fourier transforming with respect 
to the relative time t 

/+oo f + <*> 

dt j dxe ,m e-' kx G RyF (t,T,x), (2.11) 
-CO J — CO 

and 

/+<*> Ah r+oo 
~G RjF (CD,T,k)= dte ,at G R , F (t,T,x = 0). (2.12) 
-co _/t J — oo 
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We can then define a time-dependent spectral function 

p((0,T,k) = -2lmG R ((0,T,k) , 
and a time-dependent temperature 



e{T)= l -\mG F {co,T) 



f dlmG R {co,T) 



-l 



(2.13) 



(2.14) 



ffl=0 



This notion of effective temperature was also considered in [2] to assess thermalisation following 
an anisotropic deformation of the boundary. 



3. Setup of the model 

As a holographic model for a thermalising field theory, we consider a three-dimensional, thin 
shell AdS-Vaidya spacetime. For R2 > Ri > 0, it has the metric 

ds 2 = - [(r 2 - R 2 ) - 6 (v) (R 2 2 -R 2 )]dv 2 + Idvdr + r 2 dx . (3.1) 

The spacetime structure is depicted in Fig. 1 and it represents an infalling shock wave of 'null' dust 
in a black brane background that collapses to form a heavier black brane. The coordinate v is an 
Eddington-Finkelstein coordinate such that v = corresponds to the location of the shock wave. 
The coordinate r is a radial coordinate such that r = 00 is the planar boundary of the asymptotically 
AdS spacetime and r = is the black brane singularity. 



singularity (r = 0) 




Figure 1: The causal structure of the thin shell AdS-Vaidya spacetime. 

The metric (3.1) solves Einstein's equations with a negative cosmological constant and a delta 
source stress-energy tensor on the shock wave. Outside the shock wave (v > 0), the metric is given 
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by the BTZ solution with radius R2 



ds 2 = -(r 2 -Ri)dt 2 + n 



dr 1 



R 2 



+ r 2 dx 



with 



2R 2 



In 



while inside (v < 0), we recover the BTZ metric with radius R\ 

dr 2 



ds z 



-(r z -Rj)dt z + 



R 2 



+ r 2 dx 



with 



2tf, 



■In 



r-Ri 
r+R 2 



Ri 



r+Ri 



(3.2) 



(3.3) 



The dual picture of this spacetime on the boundary CFT is that of a field theory that is initially in 
equilibrium at a temperature d\ = ^ , which corresponds to the Hawking temperature of the initial 
black brane. After a sudden injection of energy at the time t = 0, there is a period of thermalisation 
until the field theory reaches a new equilibrium at a higher temperature 62 = || - 

We will now probe the equilibration dynamics of this holographic model using a geodesic 
approximation of two-point functions. While here, as in [1], we will focus on timelike two-point 
functions in order to be able to derive the spectral function, spacelike Green's functions as probes 
of thermality have been considered in [5, 6, 7]. Spacelike geodesies connecting two equal time 
boundary points separated by a distance i are conjectured to also compute the entanglement entropy 
of a spatial region of size £ in the boundary two-dimensional CFT [8, 9, 5, 6, 10]. Other non-local 
probes of thermalisation derived from the entanglement entropy, namely the mutual and tripartite 
information, were discussed in Federico Galli's talk at this workshop [11] and are studied in [12, 
13]. 



4. Geodesic approximation of two-point functions 

Consider a scalar operator €?(t,x) with conformal dimension A in the dual CFT. The time- 
ordered two-point function (T{0(t2, X2)0{t\,x\)}) is given by a path integral over all paths & 
that connect the two insertion points (t\,x\) and (^2,^2) on the boundary: 

(Tm 2 ,x 2 Mt U M)}) = f ®&e- A ^ with Jgf(£») = |^^v^^dA. (4.1) 

For large A, we can use a saddle point approximation, in which the sum over all paths can be 
approximated by a sum over all geodesies connecting the boundary endpoints [4] 

(7W2,* 2 )^i,*i)}>~ £ e-**, (4.2) 

geodesies 

where ££ denotes the geodesic length. However, due to contributions near the asymptotically 
AdS boundary, the geodesic length between two boundary points contains a universal divergence 
and needs to be renormalized. Throughout this section, we will define a renormalized length as 
8J£ = Jz? — 21nro, in terms of the bulk cut-off ro, by removing the divergent part of the geodesic 
length in pure AdS. The renormalized two point function can then be approximated by 

{T{0(t 2 ,X2)0(h,xi)}) m ~ e- ASy , (4.3) 

where 8J>f is the renormalised length of the geodesic that connects the points (?i,jci) and (?2 5 *2) 
on the boundary. 
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4.1 Continuation to Euclidean signature 

While for static spacetimes a Wick rotation to Euclidean signature is straightforward, this is not 
the case for time-dependent geometries as the AdS-Vaidya background (3.1). Therefore we perform 
a non-standard analytic continuation of the metric first proposed in [1]. Let E > R2 > Ri > and 
consider the auxiliary 'spacelike' Vaidya metric 

ds 2 = _ [(r 2 _ R 2 ) _ e(v/E)(R 2 2 -R 2 )] dv 2 + 



(r 2 -R 2 )-6 (v/E) {R 2 2 -R 2 )+E 2 
dr 2 



(r 2 -R 2 )-d(v/E)(R 2 -R 2 )+E 2 

Through the coordinate transformation 



+ r 2 dx. {A A) 



t - ^arctanh ( f y] 1 + ^ + ^arctanh (f ) for v > 
t - ^arctanh ( f \/l + ^^J + ^arctanh (§-) for v < 



(4.5) 



we recover the BTZ metric (3.2) with radius R2 for v > and BTZ metric (3.3) with radius R\ for 
v < 0. The coordinate transformations are such that v = describes the geodesic of an infalling 
shell of 'spacelike' particles in the BTZ spacetime. Observe that for E > 0, 8(v/E) = 8(v), so that 
in the E — > 00 limit (4.4) reduces to the Lorentzian 'null' Vaidya metric (3.1). 

On the Lorentzian 'spacelike' Vaidya metric (4.4), we can now perform an analytic continua- 
tion on the time coordinate w = iv, as well as on the parameter Q = iE. Without loss of generality, 
we can take Q > and find the metric 



ds 2 = [(r 2 -R 2 )-6(w/Q)(R 2 -R 2 )] dw 



2Qdwdr 



r 2_ R 2 ) _ e(yv/Q) ( R 2_ R 2 ) _ Q 2 
dr 2 



{r 2_ R 2 ) _ e{w/ Q ){R 2_ R 2 ) _Q2 



+ r 2 dx 2 , (4.6) 



where as before 0(w/Q) = d(w). Note that the radial coordinate r now runs from \j Q 2 +R\ to 00. 
Letting 

/(r)2s {/-W = ^-2 2 -*? fo-<0 
' fcm(r) 2 = r 2 -Q 1 -Rl for»->0 



the metric becomes 

d* 2 = f(A 2 dv,! 2 J- ( Odu! - 

f(r) 



dr \ 2 

ds 2 = f(r) 2 dw 2 + ( Qdw - — T ) + r 2 dx 2 , (4.8) 



which is manifestly positive. 

4.2 Geodesic length in Euclidean Vaidya 

We now compute the length of the geodesies that connect the points (ro, T\ , x\ ) and (ro, 12,^2), 
where ro denotes the location of the regularised AdS boundary. We will assume At = T2 — Ti > 
and, for computational simplicity, mostly focus on equal space geodesies with Ax = X2 — x\ = 0. 
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The geodesies in Euclidean Vaidya consist of piecewise geodesic curves in BTZ. Above the shell 
(w > 0), we have the geodesies [1] 



w±{r) 



xn ± — arctanh 
R 2 



To ± — arctan 
R 2 



A±(r) = Ao±arccosh 



and below the shell (w < 0), we have the geodesies 




(4.9) 



Ri 



-arctan 



^2 

Q 



(4.10) 



(4.11) 



x± r 



w±(r) 



xo db — arctanh 
R\ 



Tn ± — arctan 



A±(r) = Ao±arccosh 




(4.12) 



1 t (Rx 
— arctan — 

Ri \Q 



(4.13) 



(4.14) 



Here db denote two separate branches (appearing above and below the turning points of BTZ 
geodesies) and the parameters appearing in these expressions satisfy the conditions ^ T 2 ^ R\ < 
^andO^f 2 ^R 2 <t 2 + . 



r 2 and ^ r 2 ^ R 2 < T 2 . We can now distinguish three cases. 



1) If both endpoints are above the shell (T2 > Ti > 0), then there is a geodesic connecting them 
which is entirely above the shell. Any additional geodesies would require the existence of 
geodesies below the shell in BTZ that connect two points on the shell, and, as we show below, 
these do not exist. 

2) If both endpoints are below the shell (0 > T2 > Ti), then there is a geodesic connecting them 
which is entirely below the shell. Any additional geodesies would require the existence of 
geodesies above the shell in BTZ that connect two points on the shell, and again these do not 
exist. 

The geodesic length is then given by the BTZ result, derived for example in the Appendix of [1], 



<5j^ = ln 



R 2 



sm 



RAt 



+ sinh 2 



RAx 



(4.15) 



where R is respectively R 2 and Ri in case 1) and 2). 
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The existence of geodesies that have two endpoints on the shell, and that are above (or below) 
the shell would require the existence of a local maximum (or minimum) of w(r) (or vP(r)). We note 



dw 4 



other hand 



r Q ) = has no real solutions, so no geodesic in BTZ can have a local maximum. On the 

dw- (_ \ n i u., _ gr + r_ 



-5— (r P) ) = can be solved by r n - 7=—, : — 

we consider here, all have r = and thus no local minimum 



,2_ f 2 



--. Equal space geodesies as 



3) Finally we consider geodesies with one endpoint above and one below the shell (T2 > > Ti). 
They cross the shock wave once and hence we need to determine how they refract at the shell 
location in order to compute their length. 

The conditions that need to be imposed for the geodesies at the shell (w = 0) were derived in [1]. 
They are given by the following refraction law 



1 



f2dx 
J dr 



in 



Ar [ Qf 

f2dx 1 ^ J 



(/£)' 



dr 
1 



dw 
dr 



1 



(/if 



, (4.16) 



which is supplemented by the continuity conditions x- m {r*) = x ou t{ r *) an d w,- n (r*) = w out (r#) = 0, 
where is the value of the radial coordinate at which the geodesic reaches the shell. 
The boundary conditions at the AdS boundary are 



x+(r- 

w + {r - 

X-(r - 

W-(r 
v v 



' °°J = X2 

> °°) = %2 
•00) = X \ 

> 00) = % x 



X2=xq + ^arctanh 



T2 = Tq + ^arctan 



(aj-r 2 -: 



x\ = xq — ^-arctanh , , 



(4.17) 



T\ = T - 



arctan 



(r 2 + -« 2 ) 



and the continuity conditions 



1 u r 

xq + — arctanh I — , 



: -r 2 + 



1 u I r 

■ xq H arctanh 

R\ 




(4.18) 



To H arctan 



To H arctan 

R\ 




1 [R X I {(?+R\ 

+— arctan | — \j 1 



Q 



1 f Ri 
— arctan — 

R2 \Q 



1 . (Ri 
— arctan — 

Ri \Q 



(4.19) 



(4.20) 
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In terms of the solutions (4.9)-(4.14), the refraction conditions read 

r 2p2 p2p2 



R 2 R \ 



(4.21) 




R2\l{rl-Tl){rl-T\)-Qr* 



rl-Q 2 -R\ 



Rx\l{rl-T 2 _){rl-T\)-Qn " ' — 1 . (4.22) 



R\-Tl)(T\-R\) 



Altogether, these are nine conditions, for nine unknowns (xo, To, xo, To, T + , T + , T+, r and r*). 
The renormalised length of the geodesic is given by 

8J£ = hm (A+(r ) - A_(r ) -21n(r )) = Ao - Ao + ^ln U^fA +\ ln ( f2 Zp ) > 

(4.23) 

which, using the continuity condition A + (r*) = A + (r*), can be written as 



6^ = In V , v , . (4.24) 

We will focus on equal space correlators, Ax = x% — X\ =0. Because of the refraction condition, 
^ does not change sign at the shell. Also the geodesies in BTZ have no local extremum of x(r). 
Therefore, we must set F = F_ = 0. It then follows that xq = xq = x\ = x 2 . We are thus left 
with five conditions for five unknowns. After solving this system, we find that the renormalised 
geodesic length is given by 



8JZ = In I == | , (4.25) 

_ (r, - cos(7) J rl - Rj)(r* + cos( 7 ) Jrl -R\) 



where the unknown r* is implicitly determined by the relation 1 

1 I , s 2sin(r) \ 1 / <2sin(y) , 

' cos( 7 ) - . g Ul = . cos(r) - . g Ui , (4.26) 



\Jr 2 *-Q 2 -RlJ V r *-^V ^rl-Q 2 -R 2 K 
and we have used the expressions 

sin(/(r, ) ) = 2121 , (4.27) 

r m -co&(R 2 n(r m ))Jr*-F% 



'The expressions (4.26)-(4.28) were obtained assuming ^ y(r*) ^ tt/2 and ^ Y{r*) ^ Tt/2. Outside this range, 
similar expressions can be derived from (4.22). 
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and 



with 



and 



sin(/(r*)) 



Ri sin(7?iQ(r*)) 



n-cosiRyQir^Jrl-Rj 



R 2 



-arctan 



H(r*) = Ti + —arctan 
^1 




(Q 2 + R\) 



R 



-arctan 



QJ 



(4.28) 



(4.29) 



(4.30) 



The parameters T + and T + in the geodesies (4.9)-(4. 14) are given by the formulas T 2 (r*) =R\ + 
(r 2 sin 2 (y(r*)) and f^r*) =R\ + (r 2 - /? 2 ) sin 2 (y(r*)). In Fig. 2, we plot the profile of a 
sample of equal space geodesies in the Euclidean AdS-Vaidya geometry. 



w 




Figure 2: Equal space geodesies in tlie Euclidean shell background for R\ = 0.5, Ri = \, Q — 0.25, Ti = — 1 
and, increasing from the bottom up, Ti = —0.25,0.25,0.5,0.75, 1, 1.25. 

If we perform the double Wick rotation Ti = it\ , T2 = ih and Q = iE and take the limit E — > 00, 
in order to recover the 'lightlike' Vaidya result, we finally find 



8^ = In 



- . (Rit 2 \ . , fR\i\ 



cosh 



sinh i 



2 . .{R2h\ .fRih 
■ sinh 1 — — I cosh I — — 



(4.31) 



r 1 v 2 / V 2 / R 2 V 2 

which reduces to the expression obtained in [1] in the R\ = case. 

4.3 Two-point functions 

4.3.1 Equilibrium CFT 

The Euclidean two-point function associated to a CFT in thermal equilibrium at a tempera- 
ture 6 = R/ (2k) can be computed in the geodesic approximation from the renormalised geodesic 
length in Euclidean BTZ (4.15). Remember that, in the geodesic approximation, the renormalised 
Euclidean two-point function in the dual CFT equals 



T 2 ,x 2 )0(Ti,x l )) rm ~ e 



(4.32) 
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From (4.15), we then find the (Euclidean) thermal two-point function 



^thermal /_ . — 



,2 fRAr^ 



,2 (RAx 



w v snr {'-^) -l-sinh" 



)))' 



(4.33) 



Despite the fact that we have used a geodesic approximation, the result for the thermal two-point 
function is exact, being fully constrained by conformal invariance (apart from an overall scaling). 
Note that we can obtain the vacuum result by taking R — > 



^vacuum 



(T 2 ,X2;T hXl ) = lmiGl™ 1 ^,^,*!) = — — 

(Az 2 +Ax 2 ) 



(4.34) 



It is well known that by Wick rotating the Euclidean two-point function, one finds the time-ordered 
(Feynman) two-point function. Therefore let T = lim e ^o^'^ e ^ — it, such that T 2 = — t 2 + ie and 



iG F (t 2 ,x 2 ;t\,x\) = G E (it 2 ,X2\ih,xi) . 



(4.35) 



The retarded two-point function can be obtained from the time-ordered one by applying the relation 
(2.3). For non-integer scaling dimension A, 2 this leads to 



•^thermal/, v , t n 

iG F {t 2 ,x 2 ,ti,xi) = 



6 [(Ax) 2 — (A?) 2 j a — (ajcj-j e ■ — 

( S inh 2 (*f£) -sinh 2 (f^))] A + (sinh 2 (^) -sinh 2 (*fs 



6 [(Af) 2 -(Ax) 2 ] e- i7lA 



i A ' 



t0)1 

(4.36) 



and 



^thermal/, „ „ \ _ 



2 sin ttA (A?) ^— ^ K -— -J -r. (4.37) 

[£ (sinh 2 (™)- sinh 2 (^))] A 



4.3.2 Thermalising CFT 



In an analogous way, from (4.31) we can then find (the geodesic approximation of) the Feyn- 
man and retarded two-point functions in the thermalizing CFT dual to three-dimensional Vaidya: 



iGf (t 2 ,x',ti,x) 



e-^d(- tl )d(-t 2 ) 
^sinh^-?!))' 2 



+ 



e-™ A 0(h)0(t 2 ) 



^anh(^-fi)) 



2A 



+ 



+ 



e-« A d( tl )d(-t 2 ) 



2A 



r\ cosh (¥) sinh (^) - I sinh (^) cosh 
e- inA d(-h)d(t 2 ) 

2A 

I cosh {*f) sinh (^) - | sinh cosh (^) 



(4.38) 



The case of integer A is discussed in detail in [1]. 
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and 



G R (t 2 ,x;h,x) 




e(-h)e(-t 2 ) 



e(-h)d(t 2 ) 



Rl cosh(^)sinh(^)-^sinh 




In the limit R 2 — > Ri of equal initial and final temperatures, these expressions reduce to the equi- 
librium ones, (4.36) and (4.37). Also, in the limit of zero initial temperature, R\ — > 0, we recover 
the expressions for a thermalising CFT that were presented in [1]. 



5. Spectral function and effective temperature in thermalising CFT 

In a two-dimensional CFT in equilibrium at a temperature 8 = R/ (2n), we know the spectral 
function analytically [1, 14]. It is obtained from the Fourier transform of the retarded two-point 
function (4.37) using equation (2.6): 



p(C0,k) 



R 



2A-2 



/ A i(co + k) 

u + — 



2R 



A /(ft) - k) 
2 + 2R 



When averaged over momenta, it becomes 
P(0») = 



dk , 2R 2A - 1 



r a+ 



1(0 

~R 



sinh 



sinh ( 7 ^-\ 
\ R J 



K(0\ 
~R~) 



(5.1) 



(5.2) 



To compute the momentum average of the time-dependent spectral function (2.13) in a thermalis- 
ing CFT, we need to perform a numerical Fourier transform of the equal space retarded two-point 
function (4.39). As explained in section 2, one can also define a notion of time-dependent tem- 
perature (2.14) from the Fourier transforms of the thermalising Feynman and retarded equal space 
two-point functions. The results of this analysis are shown in Figures 3 and 4. 
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Figure 3: The spectral function as a function of the frequency CO for different values of the average time T. 
In both plots, we have taken A = 2.25 and R 2 = I. The plot (A) has R\ = and plot (B) has Ri = 0.5. 

In the distant past (T — > — °°) the spectral function is given exactly by the equilibrium expres- 
sion (5.2) for a temperature Q\ =R\/ (2n), while in the distant future (T — > oo) it coincides with the 
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equilibrium result for a temperature 62 = R2/(2n). The notion of time-dependent spectral function 
that we have defined interpolates smoothly between these curves. However, for a definite frequency 
(0, p(co,T) does not increase monotonically from the past limiting value to the future one, but is 
an oscillating function. 




T 



Figure 4: The temperature as a function of average time for different values of the conformal dimension 
A. In both plots, we have taken R2 = 1 and R\ = 0.5. In plot (A) we have taken different values for A as 
indicated in the legend. Plot (B) lias A = 2.25. 

For sufficiently small values of A, the time-dependent notion of temperature that we have 
defined gives a smooth transition from the initial temperature to the final temperature. However, 
for larger values of A there are values of the average time for which the spectral function will have 
zero slope at 00 = 0. This results in a divergence in the temperature just before the time T = as 
shown in the plots. Similar singular behavior of the effective temperature was also observed in [2]. 
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